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Abstract 

We prove shock formation results for the compressible Euler equa- 
tions and related systems of conservation laws in one space dimension, 
or three dimensions with spherical symmetry. We establish an L°° 
bound for C 1 solutions of the one-D Euler equations, and use this to 
improve recent shock formation results of the authors. We prove analo- 
gous shock formation results for one-D MHD with orthogonal magnetic 
field, and for compressible flow in a variable area duct, which has as a 
special case spherically symmetric three dimensional flow on the exte- 
rior of a ball. 
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1 Introduction 

One of the defining features of systems of hyperbolic conservation laws, 

Ut + f(u) x = o, (i,t)elxl+«el", (l.i) 

in the absence of regularizing dissipative effects such as viscosity and heat 
loss, is the blowup of gradients and subsequent formation of shock waves. 
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Once a shock forms, discontinuous weak solutions must be considered, and 
the corresponding loss of regularity makes the systems hard to analyze. 

We are concerned with the conditions under which a shock forms in a 
classical C 1 solution of (jl.ip . For pairs of equations (n = 2), Lax showed 
in [131 HH] that shocks always form for data which is anywhere compressive. 
In work initiated by John, this was extended to larger systems, given restric- 
tions on the size of the initial data H6| 121] . In particular, all waves in 
those solutions are assumed to be weak. One of our goals is to remove these 
severe restrictions on the initial data, which includes allowing for large waves 
in the solution. We restrict our attention to systems arising from physical 
problems because of the variety of unstable phenomena which can affect 
general 3x3 systems when the data is allowed to be large [10J [271 1211 SO] . 

Recently, the authors have obtained large data results for the 3x3 com- 
pressible Euler equations [21 E]- In all of these works, (jl.ip is differentiated 
and manipulated to get a system of PDEs with quadratic nonlinearities for 
the evolution of gradients along characteristics. For two equations, these 
PDEs decouple to yield a pair of Riccati type ODEs along characteristics; 
the quadratic nonlinearity then implies that the derivatives must blow up in 
finite time |131I14|. For general systems, if the data has small total variation, 
the nonlinear waves separate in finite time, after which they are essentially 
decoupled and genuine nonlinearity in each field causes gradients to blow up 
by the same mechanism. 

In [2113], we consider the compressible Euler equations, and make a series 
of nonlinear changes of variables to get an uncoupled pair of quadratic ODEs 
for the genuinely nonlinear fields, which exhibit the same blowup mechanism. 
Although there are no restrictions on the size of the initial data, the results 
of [2113] rely on assumed global bounds for the state variables. In this paper, 
we first remove this restriction by obtaining an a priori estimate on the 
state variables for C 1 solutions of the Euler equations, yielding a condition 
on only the initial data which guarantees blowup of gradients. See [19J for a 
description of attempts to obtain a stronger L°° bound and global existence 
which were unsuccessful. 

We then extend the method to larger systems of physical interest. We 
show similar gradient blowup results for one-dimensional magnetohydrody- 
namics (MHD) with orthogonal magnetic field, and compressible flow in 
a variable area duct. As a special case, we obtain a gradient blowup re- 
sult for spherically symmetric three-dimensional Euler flow on the exterior 
of a ball. These blowup results again rely on global bounds for the state 
variables, but in view of our obtained bounds for the compressible Euler 
equations, we expect that this is just a technical assumption. Indeed, we 
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note that all known cases of blowup of state variables occur in non-physical 
systems HHJ E7J EH ED] ■ 

The mechanism for gradient blowup and shock formation is the steepen- 
ing of compressions due to genuine nonlinear ity: if the wavespeed is larger 
behind the wave, it will eventually break. For 2x2 systems (or diagonal sys- 
tems such as chromatography), the different wave families weakly decouple 
and all compressions will collapse into shocks. However, in 3 x 3 systems, 
and especially for large data, the wave interactions between different fami- 
lies influence this process and can delay the onset of shocks p2J EUJ [25] . In 
the context of Euler, the contacts act as partial reflectors of waves and with 
care one can set up solutions in which compressions are cancelled before 
they focus [251 [5]. 

Our conditions which guarantee blowup are expressed in terms of global 
bounds for the variables, so we must ensure that the state does not leave 
some compact set. In particular, for a 7-law gas, the vacuum must be 
avoided. Moreover, when 1 < 7 < 3, our estimates depend on a lower 
bound for the density as well as on the upper bound. This is also true in 
the isentropic case treated by Lax. 

The paper is organized as follows: in Section 2, we obtain a priori L°° 
bounds for the density and velocity in the one-D compressible Euler equa- 
tions. We then use this to stregthen the gradient blowup result of [2]. In 
Section 3, we consider one-D MHD with orthogonal magnetic field, derive 
Riccati type equations for gradients, and use these to obtain a stronger sin- 
gularity formation result than that of [4]. In Section 4, we similarly derive 
Riccati type equations and present a shock formation theorem for com- 
pressible Euler flow in a variable area duct. This includes the special case 
of three-D spherically symmetric flow in the exterior of a ball. 

2 Compressible Euler equations 

We begin by considering the initial value problem for the compressible Euler 
equations in a Lagrangian frame in one space dimension, 

74-^ = 0, (2.1) 

u t + Px = 0, (2.2) 
^u 2 + e) t + (up) x = 0, (2.3) 

where p is the density, r = p~ x is the specific volume, p is the pressure, u is 
the velocity, e is the internal energy, and x is the material coordinate. The 
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system is closed by the Second Law of Thermodynamics, 

TdS = de+pdr, (2.4) 

where S is the entropy and T the temperature. For C 1 solutions, it follows 
that (|2.3p is equivalent to the "entropy equation" 

St = 0. (2.5) 



When the entropy is constant, the flow is isentropic, and (|2.ip and (|2.2j) 
become a closed system, known as the p-system |24| . 
We assume the gas is ideal polytropic, so that 

pr = RT and e = c T T 

with ideal gas constant R, and specific heat c T ; this implies 

p = Ke^T~~< (2.6) 

with adiabatic gas constant 7 = — + 1, see [fjj. The nonlinear Lagrangian 
sound speed is 

c : = V-p r = a/^t" 2 ^" . (2.7) 

2.1 Coordinates 

We use the coordinates introduced by Temple and Young in [25] • Define 
new variables m and r] for S and r, by 

m — d > 0, (2.8) 



and, referring to ([27 



It follows that 



V := I ^dr= 2 -^r- 2 ^ >0. (2.9) 

m ' 



2 

t = K T rj 1- 1 , 

9 -22- 

p = K p m 2 rj-y- 1 , (2.10) 
7+1 

c = c(ry, to) = .fr c mrj-i- 1 , 
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where K T , K p and K c are positive constants given by 

K T = (^)^ : K p = KKr\ and K B = y/Kj K'^ , (2.11) 
so that also 

K p = ^-K c and K T K C = 2=1. (2.12) 

In these coordinates, for C 1 solutions, equations (|2.ip - (|2.3p are equivalent 
to 



i]t + -u x = 0, (2.13) 
m 
P 

ut + mcr] x + 2 — m x = 0, (2-14) 
m 

m t = 0, (2.15) 

the last equation being (|2.5p . which is equivalent to (|2.3p . Note that, while 
the solution remains C 1 , m = m(x) is given by the initial data and can be 
regarded as stationary. 

2.2 L°° bounds for C 1 solutions 

Our first goal is to provide a uniform upper bound for the density p and 
velocity \u\, depending only on the initial data, for C 1 solutions of (|2.ip 
(|2.3p . Denote the Riemann invariants by 

r:=u — mr], s:=u + mrj. (2-16) 

For isentropic C 1 solutions (m constant) the Riemann invariants s and r are 
constants along forward and backward characteristics, respectively. How- 
ever, for general non-isentropic flow, s and r vary along characteristics. 
The forward and backward characteristics are described by 

dx dx 

— = c and — = — c, 

dt dt 

and we denote the corresponding directional derivatives along these by 
d+:=f + c|: and a_:=&-c£, 

respectively. Using ([236]) , (|2~T0j) and fl2~T2|) . equations (I2TT31) and (I2TT41) 

become 

d + s = ^^(s-r), (2.17) 

j f) rn 

d-r = - — (r-s), (2.18) 
27 m 
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while by (|2.3[) . we have 

<9+m = cm (x), and d-m = —cm'(x). (2-19) 
Applying an integrating factor and setting 

s = m 2 t s and r = m 2 t r, 

we obtain the system 

8+?= _J_^ ? , (2 .20) 
27 m 

„ ~ 1 cLm „ 

d_r = s. 2.21 

27 m v ; 

We obtain upper bounds for C 1 solutions of f)2. 1 j) — f)2.3j) by bounding \s\ 
and \r\. We assume that the initial entropy S(x) is C 1 and has finite total 
variation, so that 

1 r +co \m'(x)\ 

V.= — \S'(x)\ dx = / illLfli dx < 00, (2.22) 

2c T y.oo y_oo m(x) 

while also, by (|2.8p . 

< M L < m(-) < M v , (2.23) 

for some constants Ml and M;y. Also, since p and |w| are bounded initially, 
there exist positive constants M s and M r , such that, in the initial data, 

|?o(-)| < M s and |r (-)| < M r . (2.24) 
We define two useful constants by 

JVi := M s + V M r + V (V M s + V 2 M r ) e y2 , 
N 2 := M r + V M s + V (VMr + V 2 M s ) ^ , 

where V := 5-, which clearly depend only on the initial data. 

Theorem 2.1. For C 1 solutions of [2. 1\) - I12.3\) . we have the a priori bounds 

, . N 1 +N 2 _L Ni + N 2 -L-1 

\u\ < — — — -Mu^i and p < — — — -M L ^< , 

Note that the density p need not be bounded away from zero, as vacuums 
can form in infinite time [26\ [5] . The theorem is an immediate consequence 
of the following estimate. 



6 



Lemma 2.2. For a given point (xi,ti), suppose the solution is C 1 in the 
characteristic triangle bounded by the forward and backward characteristics 
through {xi,t\) and the line t = 0. Then 

\s(xi,t{)\ < Ni and \r(xi,ti)\ < N 2 . 

Proof. We prove the bound on the bound of r(x\,ti) is obtained 

similarly. Referring to Figure [TJ we denote the forward and backward char- 
acteristics through a point (z* , t* ) by 

~^x* = {{ x ^i x )) I x — x *} = { (lz^(t) , t) | t < i*} and 
Ac. = {( x i*t( x )) I x > x *} = {(^(*)>*) I t < 
respectively, parameterized by x or t, as convenient. 




^i(O) ^(0) ^(0) ^i(O) 



Figure 1: Characteristic triangle 

We integrate (|2.20p along the forward characteristic ~~& Xl from (a^i(0),0) 
to (xi,t\), with a change of integration variable, and use x a as the new 
parameter. We then integrate (|2.2ip along the backward characteristics 
C Xa1 from to t a , for each {x a ,t a ) G ~C Xl , and again change variables. We 
have 

s(xx,ti) = s (^i(0)) - — r — {x a )r(x a , t(x a )) dx a , 

27 J m 

^i(O) 

where we have used (|2.19p . and similarly 



r(x a ,t a ) = roOF CT (0)) - — 7, — (x) s(x, i CT (x)) fix. 

27 y m 

^r(O) 
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Here J^and ^indicate the direction of the characteristic along which the 
integration is performed. Combining these, we get 

a(3i,ti)=3bC2i(0))-^- k — Mr fa(P)) dx a (2.25) 

27 J m 



^i(O) 



■t-i 



4 7 2 



^i(O) 



m 



+ -r-~ T — (x a ) ( kL — (x)s(x, t a (x))dx ) dx a 



in 



V CT (0) 



The first two terms can be estimated by our initial bounds, and we will 
apply a Gronwall inequality to estimate the third term. 

Equation (|2.25|) continues to hold for any point (x£,ig) on the backward 
characteristic £ Xl , with x\ replaced by x^. We multiply (|2.25j) by ^p(z^), 
take absolute values, integrate in time along the backward characteristic 
C Xl , and change to spatial variables. Doing this, we get 



Fi(0) 



Vi(0) 



| s o(^(0))| dn{xt 



Xl 

^i(O) x £ 

|f (lFa(0))| d^l(Xa) d(J,(X£) 

x i ~^(0) 
^i(O) x s ^V(O) 



1 



(2.26) 



^1 "^(0) No- 
where ci/i is the Stieltjes measure 

d/i(x) :- 



m 



m 



dx, 



which has total mass at most V along any characteristic. Again the first 
two terms of (I2.26P can be estimated by the initial data, so we focus on the 
last term. 

For {x a ,t a ) on the forward characteristic ~C Xl , set 



^t(O) 



F(x a ) := ^( \s(x, dfi(x), 



(2.27) 



S 



where the integral is along the backward characteristic C XcT . If the back- 
ward characteristic through (x, *¥a{x)) is extended up to ~C X1 , the innermost 
integral in (I2,26h is estimated by 

\s(x, ta(x)) \ dfi(x) < F(x a ), 

provided (x-&, t w ) S £ XfT , so that also tFa(O) = *x a (0), and we change the 
order of the other two integrations. Using (I2.27j) with (12.221) and (|2.24j) . it 
follows that (|236|) yields 

"5a (0) 

F( Xl ) <VM s + ^-V 2 M r + ^ Jj* F(x a ) dfx(x) dn{x a ) 

^i(O) x <? 

< V M a + ^V 2 M r + ly F(x a ) dfi(x a ). 

^i(O) 

We now use Gronwall's inequality to get 

v 2 

F( Xl ) < (VM s + ±V 2 M r )e^ =:Nt. 

It follows similarly that for (x a ,t a ) 6 C Xl , we have F(x a ) < JVi. Finally, 
using this estimate in (|2.25p , we get 

\s(x 1 ,t 1 )\<M s + ±VM r + 1 ± T VN 1 = N u 

and the proof is complete. □ 

2.3 Singularity formation 

In [2], we introduce gradient variables 

3(3-7) 7+1 „ 

y = m 2(37-1) ^2(7-1) ^ u _|_ m _ _____ mx yy) anc [ 

3(3-7) 7+1 „ 

q = m 2(37-1) ^2(7-1) ((■_ - mrj) x + g^m^), (2.28) 
and derive Riccati type equations for their evolution. 
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Lemma 2.3. [2] For C 2 solutions of \2. l \) - \2.3\) , we have 



d+y = a + a 2 y 2 , 
d-q = a + a 2 q 2 , 



(2.29) 



where 




91 3(3-7) 3(7+1) 

m x \ m 2 (37-i) rjWTil 



+1 



3(3-7) 7+1 



(2.30) 



a 2 = -K c m 2 ^- 1 ) r/ 2 ^- 1 ) < 0. 



Furthermore, 



y\ or \q\ oo z# 



or |r x | — )• oo. 



(2.31) 



Similar equations are derived independently in [18]. Using these equa- 
tions, a singularity formation result was proved in [2]. By applying Lemma 
12.21 we can improve that result to obtain our main theorem on the Euler 
equations. 

Theorem 2.4. Assume the initial data are C 2 with entropy having bounded 
variation, and suppose there is a positive constant M* such that the initial 
entropy satisfies \m"(x)\ < M*. // 1 < 7 < 3, we also assume that the 
density has a positive global lower bound. There exists positive constant N 
depending only on the initial data, such that, if the initial data satisfies 



then \u x \ and/or \t x \ blow up in finite time. 

This result is an extension of Lax's singularity formation result in [13] for 
2x2 strictly hyperbolic systems. We recover the isentropic case by taking m 
constant, so ao = 0. The conditions which guarantee blowup are expressed 
in terms of global bounds for the variables, so we must ensure that the state 
does not leave some compact set, and, in particular, the vacuum must be 
avoided. When 1 < 7 < 3, our estimates depend on a lower bound for the 
density as well as on the upper bound. This is also true in the isentropic 
case treated by Lax. However, when 1 < 7 < 3, we expect a more refined 
analysis to show that even if the solution approaches vacuum, the gradient 
still blows up in finite time [29 1. 

Proof. By Theorem 12.11 and (|2.9p . we know n has a global upper bound 
depending only on the initial data, denoted E\j. 




(2.32) 



10 



By (|2.30p . it is easy to calculate that, if > 0, 
which implies the uniform bound \f — ao/a 2 < ^V", where 



N :-- 



^g^M.Ep^Mp, l< 7 <5/3, (2 ^ 3) 



37-1 37-5 
37-1 



l{l+ l^-i) M * E U M v \ 7>5/3 



It follows that if y < —N, then ao + a 2 y 2 < 0. 

Now suppose that (|2.32|) holds. Then there exist e > and xq such that 

y(x ,0) < -(l + e)iV, 

say. Now consider the forward characteristic starting at (xo , 0) . By (|2,29D , 
along this characteristic (parametrized by t) we have d+y < 0, so also 

y(t)<-(l + e)N for t > 0, 

which in turn implies 



l/(*) 2 
(1+0 

for all t > 0. Now (12T2H implies 

<9+y = a + a 2 y 2 < (1 - jj^yz) a 2 y 2 < 0, 
since a 2 < 0. Integrating, we get 

t 

w^w>-f (1 -^ )a2dt ' {2M) 



where the integral is along the forward characteristic. By (|2.3(jp . when 
7 > 3, a 2 is negative and bounded above, so the right hand side of (|2.34p 
approaches zero in finite time. This implies that y(t) approaches —00 in 
finite time, so that \t x \ and/or \u x \ blow up. 

When 1 < 7 < 3, and assuming 77 has a positive (global) lower bound 
along the characteristic, then a 2 is again negative and bounded above, so 
the result follows similarly. It is clear that the same argument holds along 
a backward characteristic if inf q < —N. □ 
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3 Magnetohydrodynamics 

3.1 Prior results for generalized Euler equations 

We recall the results of [3j for the generalized p-system (which includes the 
smooth Euler equations) in one space dimension, 

n ~ u x = 0, 
u t + P{r,x) x = 0, (3.1) 

where P(t, x) is a C 3 function of r > and x, satisfying 

P T < and P TT > 0, 

so that the system is hyperbolic with wavespeed 

c := c(t,x) = \J—P T . 

It is convenient to introduce new variables (h, fi) for (r, x), by setting 

/T* AT* 
c dr = / \J —P T cLt and /U = x, 

where t* is a constant or infinity. It follows that 

h T = -c, Th = , and = c, (3.2) 

c 

and furthermore, for any C 1 function f(r,x), 

fh = ~— and f fl = —f T + f x . (3.3) 

c c 

We define gradient variables by 

p 

y ■= Vc (u + h) x H — p — /, 
Vc 

<?:= yfi( u -h) x -^+I, (3.4) 
Vc 

where 

I = I(h,n) = [ \V^(—) dh, (3.5) 

J h 2 \cy h 

and /io is a constant. 
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Lemma 3.1. For C 2 solutions of \3. 1)) . we have 

d + y = a + aiy + a 2 y 2 , 

cLg = a - aiq + a 2 q 2 , (3.6) 

where 




For arbitrarily given positive constants Ai and Bi, we denote by fC the 
compact set whose interior tC° is given by 

|A|<-Bi, A 2 <c<B 2 , A 3 <c h <B 3 , \ Cfl \ < B 4 
|c MM | < B 5 , \c hfl \ < B 6 , \p^\ < B 7 , |p w | < B 8 . 

Lemma 3.2. OV There exists a constant N > depending only on fC, such 
that, if the C 2 initial data of 113.1)) satisfy 

(u(x,0), r(x,0)) G /C°, /or a// x, 

and i/ 

inf { yO.O), q(-,0) } < -N, 
then there exists T* = T»QC,N) such that either 

max{|uj;|, \t x \} — > oo as t — > T*, 

or i/iere is some poini (x b ,t b ) wit/i i& < T* swc/i i/iai 

(u(x b ,t b ), T(x b ,t b )) G 3/C. 

3.2 Refinements for MHD 

In this section, we apply and upgrade Lemma 13.21 for one-D Magnetohydro- 
dynamics. The system models the motion of a smooth compressible fluid 
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coupled to a magnetic field H = (H\, H2, H3). In a Lagrangian frame in one 
space dimension, C 1 solutions satisfy the quasilinear system 



dr du\ 
dt dx 



0, 



at ox ax 



"" L +|-(p + J + *!))=<>, (3i 



at dx 



du-j dH; 



as 

at 



0, 



Here the fluid quantities are the density p, specific volume r = p , veloc- 
ity field (ui,U2,us), entropy S and hydrostatic pressure p = p(r,S), while 
(Hi, H2, H3) is the magnetic field and po is the magnetic constant [T7J. In 
one-dimensional MHD, the first component H\ of the magnetic field is nec- 
essarily constant. 

System (|3.8p simplifies significantly if we take the constant H\ to vanish, 
so the magnetic field is perpendicular to the direction of motion. In this 
case, the full system (j3.8|) can be written 

dr du\ 

dt dx 
du\ dp 

regarded as a generalized p-system, and coupled with 
dH j duj dS 

where we have written 



, 1 H$ + H? 

Hj := t Hj and p := p + - p ^ 

2 T 2 



3 



System (|3.9p is clearly an example of (|3.ip . so all results in [4] apply 
directly. Here we obtain a stronger singularity formation result than Lemma 
by restricting our consideration to a polytropic ideal gas, 



P 

with adiabatic gas constant 1 < 7 < 2. 
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Theorem 3.3. For system \3.9ji . \3. 10\) , assume that the initial data is C 2 
with uniform bounds 

|G(-,0)| < N , \G(-,0) X \ < Ni, and \G(-,0) XX \ < N 2 

for each ofG = S, H 2 , H%. Assume also that the density is globally bounded, 

0<N pL < p(x, t) < N pU , for any (x, f)elx R+ (3.11) 

for constants N p L and N p jj. There exists a positive constant N , depending 
on Nj and N p u, but not on N p L, such that, if 

inf { y{-,0), q(-,Q) } < -N, 

then u x and/or t x blow up in finite time. Here y and q are defined in \3.J^ , 
with P = p. 

Although we assume the global bounds (|3,lip . iV is independent of the 
lower bound N p l of density. Theorem 12.11 for the compressible Euler equa- 
tions indicates that it is reasonable to expect an a priori upper bound on 
the density, and such a bound would imply that iV depends only on the 
initial data. 

Proof. As in the proof of Theorem 12. 4[ from ()3.6f) , it suffices to find uniform 
bounds for the roots of the quadratic equations 



00 + a i V + 0-2 y 2 = and ao — 01 q + a 2 < 



0. 



whose coefficients are given in (|3.7p with P replaced by p. That is, we must 
choose N so that 



2 o 2 



+ 



0,1 



2a 2 



\a 2 \ 



(3.12) 



when a\ — 4aoa 2 > 0, where we calculate 



a 1 



2a 2 



Ch 



c Vc (Pjtt +I 



and 



0-1 



2a 2 



+ 



\ a 2\ 



CVc ( Pjj, 

Ch V c 



cyfc c fPfi . 

Ch c h \ c 1 
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It follows that we can choose TV provided we find upper bounds (independent 
of N p l) for each of the quantities 







Pn 
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Ch 



, \I\, and |IJ; 



(3.13) 



we treat each of these terms separately. 

In order to simplify the calculation, we denote 

A(x):=Ke^, B(x) := ~ Mo (#f + H$), 

so that A(x) and B(x) are positive bounded functions, and A(x) is bounded 
from below by a positive constant. With this notation, 



the wave speed is 



p(r, x) = A{x) t-"< + B(x) t- 2 , 



\T I p~ T = ^-fAr-y- 1 +2Bt~ 



(3.14) 



and 



/oo /*oo i 

cdr = \/~f At-"/- 1 + 2BH dr. 

It follows by (13. lip that p, c and /i are bounded above and below. We write 

< N hL < h< N hU , 

where NhU does not depend on N p i. 

1. We first consider (^jr)h- By (|3.2p we have p^h = and by (|3.3[) . 



Pn 



C ( P; 



2 VPr 



and thus 



1 



2 W T 



PxT Pt Px PtT 



2p> T 



Now, since 1 < 7 < 2, we have 

\p T \ = 7 At^ 1 + 2Bt" 3 > Tir^ 1 , 
and since ^4 is bounded from below, {^f)h ls bounded provided 

Fi = (PxtPt ~ PxPtt) t 2j+2 
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is bounded. It is routine to compute 

F 1 = ( 1 A x + 2B x t^ 2 ){ 1 A + 2Bt^- 2 )-{A x + B x t^ 2 ){ 1 { 1 + 1)A + QBt^- 2 ), 

which is bounded independent of N p i for 7 < 2. 
2. By Q, 



Px 



Px h Pt _ 



and this is bounded by a constant independent of N p l, since 

A x t-^ + B x t~ 2 



p., 



and 



hn 



Vc (7AT-T- 1 + 2 j Bt" 3 )3 



. 7+1 _ „ 7-5 

>yA x T 2 +2B x t 2 



, , 37-1 
0(1) r — ~ 



2y/<yA + 2B tT- 2 



cir = 0(1) r 



. 7-1 



37-1 



so y/ch x = 0(1) r 4~ , where O(l) is a bound depending on x and the 
lower bound for r, and thus is independent of N p i. 
3. We similarly calculate 



C\/C 



Ca/C 



7 

2c5 



2(~ f AT-y- 1 + 2Bt-' 3 )^ 



, , 37-1 
0(l)r — ~ 



c h c T r h 2cc T 7(7 + 1)At~~ / -2 -\- QBt~ 4 
which is again bounded independent of N p i. 

4. Since the upper bounds of c, h and |(^r)/t| don't depend on N p l, the 
absolute value of I is bounded above by a constant independent of N p l. 

5. By (J33D, d33]), 



fin 



1 



h 
ho 



C H (Pn 
4y/c V c 



+ 



dh, 



so it suffices to bound and \/c( ^ )v By (|3-3 



< 



(C 2 ) T ft* 



2c 2 VH 



2c^c" 



7+1 



< 0(l)r" — 



as above. Finally, estimating as above, we obtain 

'"A =0(l)r-^, 

c J h/j, 
as required. 

Since all terms in (|3.13p are uniformly bounded, it follows that we can 
choose N so that (|3.12p holds. Moreover, N is independent of N p i. The 
remainder of the proof follows that of Theorem 12.41 □ 
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4 Compressible flow in a variable area duct 
4.1 Equations and coordinates 

In this section, we first consider inviscid compressible flow in a duct of 
varying cross section a{x'). In Eulerian (spatial) coordinates (x',t'), this 
flow is described by 

a t ' = 0, 
{ap)f + {apu) x > = 0, 
(apu) t > + {apu 2 ) x i + ap x > = 0, (4.1) 
(a p E) t > + (a p Eu + apu) x i = 0, 

a = a{x') is the cross-sectional area of the duct and E = e+ % u 2 is the total 
energy (SJ [71 [T5]. The others are the standard thermodynamic variables, 
subject to (|2.4p as usual. The first equation describes the fixed duct, while 
the others represent conservation of mass, momentum and energy, respec- 
tively. For C 1 solutions, the energy equation is equivalent to an entropy 
equation, 

+ u S x i = 0. 

We use Lagrangian (material) coordinates (x,t), defined by 

a(x ) p{x ) dx' , t = t , 
so that 

dx = a p dx' — apu dt , and dt = dt' . (4-2) 
We define a specific length by 

v := -, (4.3) 
ap 

and rewrite system f)4. 1[) in the Lagrangian frame, to get 

v t -u x = 0, 
u t + ap x = 0, (4.4) 

s t = o, 

valid for C 1 solutions, which is a variant of the generalized p-system (|3.ip . 
We restrict our attention to a polytropic ideal gas (I2.6p . so that 

p = p(v,a,S) 
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and define the generalized Lagrangian wavespeed by 

I 7-1 7 + 1 S 

C = \J—a°p v = \JK^a 2 v 2 e 2c r . 
As above, we define new variables m and z for S and v by 

7-1 



s 



m = e 2c r and z = / dv = v 2 . 

Jv ™> 7 - 1 

It follows that 

2_ „ _2t_ 

v = K T z t- 1 , p = K p a 7 m z^ 1 , 

_ 7-1 7+1 

and C = C(m, z,a) = K c a "mzi- 1 , (4-5) 

where the constants -fT are given by (|2,lip . In these coordinates, C 1 solutions 
of (|4.4p satisfy the equivalent system 

C 7=1 

2t H a 2 Ux = 0, 

m 

7_1 

ut + mCa 2~ ^ + 2 — m x — ^ pa x = 0, (4.6) 

m 

mt = 0. 

4.2 Coupled Riccati equations 

Following [2], we define gradient variables by 

7-1 ~_i 7-1 

a = + a 2 7Ti -+- J — a 2 777^, 2, 

■y— 1 -1 y— 1 

ft = u x — a 2~ 777, ^ — 2— — a 5~ 777,3, z, (4-7) 
and denote derivatives along forward and backward characteristics by 
d + = d t + C d x and cL = d t — C d x , 

respectively. 

Although the area a(x') of the duct is stationary in spatial coordinates, 
it is moving in material coordinates. We write 

da(x') , .. d 2 a(x') 
a := — and a •— 



dx' dx' 2 
so that, according to (|4.2p . a(x,t) satisfies 

at = u(x,t) d(x'), a x = v(x, t) d(x') and a xx = i> 2 a + u x a, (4.8) 

for x' = x'(x, t). 
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Lemma 4.1. For C 1 solutions of |^.i[ ), we have 

d+a = ki (k 2 (3a + (3) + (a (3 - a 2 )) + k£ (a - (3) + F(x, t), 

d-(3 = h (-k 2 (a + 3/3) + (a /3 - /3 2 )) + k^ (f3 - a) + F(x, t), (4.9) 

,± 7-1 -1 ■ j_ 3(7-l) 2 -Z±l ■ 

fc 3 = — a ± s m z a a ' 

and where 

F(x, t) = ^ 7 8 ~ftP m 2 z t- 1 a~ 7-1 (a a — 7 d 2 ) + ?n m x z 2 a -7 d. 

We note that the coefficients kj and zero-order terms F depend only 
on the solution itself and on stationary derivatives of initial entropy and 
duct shape. It is reasonable to assume that |d|, |d| and \m x \ are uniformly 
bounded, so the coefficients and inhomogeneity are essentially zero order 
terms. As in all of the simpler examples, the quadratic nonlinearity drives 
the growth of a and/or j3 to infinity, at which time derivatives blow up and 
shocks form in the solution. 

Proof. We derive the equation for a. Differentiating (|4.7p . we get 

d + a = [u xt + C (a _2 2~ mz x ) x ] + [(a _2 2~ mz x ) t + C u xx ] 
+ 2=i [(a"^ m x z) t + C (a _: V m x z) x ] . 
Now differentiate (|4.6p in x to get 

-y — 1 "f — 1 *f — 1 

(a 2 m Zx ) t + C u xx = -C x u x - (a ~ m) x z t + (a ~) t m z x , 

u xt + C (ma 2- Zx ) x = -(2 — m x ) x + ( r ypa x ) x - C x ma ~z x , 

m 

and substitute in to get 

a p 7— 1 

d + a = -(2 — m x ) x + (ipa x ) x - C x ma ~z x 
m 

-C x u x - (a _2 2- m) x z t + (a~ 1 i-) t mz x (4.10) 
+ ^ m x z) t + C (a _2 2- m x z) x ] . 

Note that we have eliminated all second derivatives except for a xx and m xx , 
since m x t = 0, and we use (|4.8p to eliminate a t , a x and a^. The quadratic 
terms in (|4.9p are those involving C x . 
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Finally, we express the derivatives u x , zt and z x in terms of a and j3. By 
([HD , dUSD and we have 

m x = ~ (a + /3), 

Substituting these into (|4,10p and simplifying yields (|4.9p . We omit the 
details, which are similar to those appearing in [2]. □ 

4.3 Decoupled Riccati equations 

As in OH], equations (|4.9p can be rewritten as a decoupled system with 
varying coefficients. Since the a equation has terms linear in j3, these can 
be eliminated by an integrating factor. By (|4.5p . (|4.6j) and (|4.7p . 

7+1 i 7-1 

d + z = z t + C z x = -KcZi- 1 ((3 + - L ^-m x za 2 ), 

so that 

1 _2±I _ y i 

p = — — z t- 1 a+z — '—m x za 2 . 

We substitute for f3 in (|4.9p . move terms including d + z to the left hand side, 

7+1 

and multiply by z 2(7-1) , to get 

7+1 , 7-1 7+1 j.+ 7+1 

9 + ( z 2(7-i) a ) + A- mx .a — 2- Z 2 ( 7-i) a +z _ WTT)d + z 

= fei(fe(3a — : ^-m a; za~ 2 2 _ ) + (— :) ^-m x 2;a~ 2 2 _ a — a 2 ))z 2 <T- 1 ) 

+ k£(a + ^m x za~^)z 1 ^) + F(x, t)^ 2 ^ 1 ) . (4.11) 

Lemma 4.2. There are variables Y and Q, defined by 

7+1 ~ 
Y = z 2 (7~i) a 4- u, m, m x ,a, a), 

7+1 ~ 
Q = z 2(7_1) /3 + Q(z, u, m, m x ,a, a), 

which satisfy the equations 

d + Y = d + d 1 Y + d 2 Y 2 , 

d-Q = d -d 1 Q + d 2 Q 2 , (4.12) 
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with coefficients 

di, di = di, di(z,u,m,m x , a, a), and 
d , d = d , do(z,u,m,m x ,m xx ,a,d,d). 

Proof. For convenience, we take 7 7^ 5/3 or 3. Writing (|4.1ip as a derivative 
along the characteristic, we get 

7 + 1 . 'V — 1 7-1 37-1 

Q (^a(-y-i) a ) + Q (—J. -m x a 2 £2(7-11) 

7(37 - 1) 

3(7 - l) 2 7 + 1 7-3 . ( 7 _1) 2 7-3 , 

v ' ' m a a a~ z 2(7-1) +z + # + / M ^2(7-1) ) 

8A C 2a c (7 - 3) 

(7 — l) 2 1 



2K C (7 - 3) 



7-1 _7^is _32_i_ (7-1) 



2 



7-3 



7 — 1 7—1 '-y — 1 7 — 1 „ . 7+1 

+ fei(/c2(3a m x za 2 ) + ( m x za a — a ))z 2 ^-^ 

7 7 

, 'V — 1 7-1 7 + 1 7 + 1 

+ k^(a + - m x za ~ )z 2 (7-i) + F(x,t)z 2 ^-^ . 

Next, we eliminate the d+u term, which comes from differentiating fcg . 
By (I4.6p . we have 

<9 + u = ti t + Cu x = —ma~ :L 2~d + z — -^-C z m x a~ x ^r + r ypa x , 

and we again write this as a total derivative plus zero-order terms in z and 
u. 

Finally, we define 

7+1 l%_1~! 2 1 7-3 I 7-1 37-1 

y = z 2(7-i) a + _^il_ yUa -iaz2(7-D + -^-l^m xa -— z 2(7-1) 
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and simplify, to get 

0+X--ZK1 ) (} + Uk c ( 7 -3)(3 7 -5) maa 2 ) 

+ 2K~fr-3) a^az 2 ^ {-^j-C zm x aT^ + jpa x ) 

-. 7 — 1 -i "V— 1 T~t~^ 

+ A:i(/c 2 (3a - 1 ^-m x za a~~) + (- :L ^m x za~^~ a - a 2 ))^ 2 ^- 1 ) 

, _-| 7 — I 7 + 1 7+ 1 

+ kj(a + ?—±m x za —)z 2 ^~ 1 ) + F(x, t)z 2 ^-^ . 
Using (|4.13p to write 

7 + 1 _ 7 + 1 _ 7+1 

Y = z 2 (->~v a + y , so a = z WF^ Y - z jjq, 
we eliminate a to get 

d+Y = do + dxY + d 2 Y 2 , 
with coefficients di containing no derivatives of z, u, and in particular, 

d 2 = -^K c z^. 

2 

Using zt- 1 = K T a p, after a tedious calculation we get, for 7 7^ 5/3 or 3, 

d = a 2 ^ (lK c K7^(^mm xx - ^±g§^ mg) a 2 p^ 1 

37-1 . 57-3 

+ L\m x au p 4 +L2mm x dp 4 



2 o 37—5 37—5 

+ L^ma a up* + L^mda up 4 

t 2—2-2 5 T~ 7 r .. — 1 2 7 ~ 3 

+ -L5 ma a p 4 +Lgaa up* 

T —2-22 7 ~ 3 r 2 •• — 1 5 T~ 7 \ 

+ L-j a aupt+L%maa p*J, 



di = Lgm x ap 2 +L 10 ua a + Lnma hp 2 , 
d 2 = -2±iK c i^ fl V p V < 0, 



where the are constants depending only on 7. 
Similarly, if Q is defined by 

7+1 i\2 , 7-3 , 7 -l 37-1 

Q = z2(7-d a+ ^ M^a^ _ _Jzl_ m;!;a -— z2(7 -d 

(3 7 -13)( 7 -l) 3 ■ -2±i *Z=^ 
+ 4K c ( 7 -3)(3 7 -5) maa Z ' 
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then it satisfies 

d-Q = d + d 1 Q + d 2 Q 2 , 

with, for 7 7^ 5/3 or 3, 

do = a*? [\K c K^{^rnrn xx - %^ ml) a 2 p^ 

37-1 - 57-3 
+ Lim x au p 4 + L2inm x a p 4 

— Q 9 37—5 — -I 37—5 

+ L^ma a up 4 + L^maa up * 

— 9 _9 2 ^7 — 7 — * Q 7 — 3 

+ L5771 a a p * + L^aa up* 

— _q 9 7—3 — q 1 57—7 , 

+ L7a aupi+Lgmaa pi), 

— — 7+1 — 1 — _-i 7— 1 

di = L$m x ap 2 +Lioua d + Lnma dp 2 , 

and are constants depending only on 7. It is clear that similar equations 
hold for 7 = 5/3 or 3. □ 

4.4 Singularity formation 

Having obtained Riccati equations (|4.12p for derivative variables, we obtain 
a gradient blowup result as in Theorem 12.41 The result states that a strong 
enough compressive wave will eventually form a shock, provided the state 
variables remain finite. Theorem 12.11 indicates that state variables should 
indeed remain bounded for all times, although a vacuum can form in infinite 
time [26|l5]. 

Corollary 4.3. Assume there exist some positive constants Mi, such that 
the duct and initial data satisfy 

Mi < a(x,0) < M 2 , |d(x,0)| < M 3 , |d(x,0)| < M 4 , 
Af 5 < m(x,0) < M 6 , \m x (x,0)\ < M 7 , \m xx (x, 0)| < M 8 , 

while the solutions satisfy the a priori bounds 

\u(x,t)\<M 9 , p(x,t)<M w , 

and p(x,t) > Mn if 1 < 7 < 3, for (x,t) E 1R x R+. T/ien there exists a 
positive constant N, depending on 7 and Mi, such that, if 

Y < —N or Q <-N 

somewhere in the initial data, then a singularity forms in finite time. 
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Proof. The roots of (|4.12p . if they exist, are given by the quadratic formula 

-di ± Jd'i -Ad d 2 di±Jdj-4:d d2 

— or 

2d 2 2d 2 

which have bounds depending on the ratios |di|/|d2|, l^ol/l^l, I di 1/1^2 1 and 
| do I /I d 2 1. It follows that these roots have a lower bound — N depending on 
the Mi, and the proof follows in the same way as Theorem 12.41 □ 

4.5 Spherically symmetric Euler equations 

The equations for compressible flow in a duct include spherically symmetric 
three-dimensional Euler flow as a special case, on the exterior of a ball. 
Using r = \x'\ as the radial coordinate, the Euler equations reduce to the 
duct equations with x replaced by r and a(r) = r 2 . 

In order to get cleaner equations, we consider the isentropic case, with 
m = 1, say, and 7 7^ 5/3 or 3. With these assumptions, m x = 0, d = 2r and 
a = 2, and (|4.12j) hold with coefficients 

r 2 ^- 2 p 2 ^ (G x u p 2 ^ 1 + G 2 p^ 1 + G 3 u 2 ) 
r- l {G A u + G^p^) 

for constants Gj. Thus the ratios are 

d\ 7-5 7-3 . ~ * 7-1 . 
— = r 2 p 4 G411IG5/) 2 
d 2 

<h = r7 _ 5 ^ {di up ^ + Q 2 pr-i + G 3 u 2 ), 
d 2 

for constants Gj. It follows that the roots (3± = dl± V^i~ 4dod2 Q f fjzt. 12|) 
satisfy the bound 

P± > -r^ p^ (G 6 \u\ +G 7 pV), 

for constants Gj depending only on the equation of state. 

A gradient blowup result similar to Theorem 12.41 now follows exactly as 
before for the exterior problem when 1 < 7 < 5. 



do = 
di = 

d 2 = 
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